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Abstract 

We study the radiative leptonic D meson decays of Dfs — > / + ^7 (I = e, /i,r), D° — ► 
uuj and D° — > / + /~7 (I = e, fi) within the light front quark model. In the standard model, 
we find that the decay branching ratios of D~)\ — > e + z/ e 7, Dt* — > y^v^ and Df\ — > r + i/V7 
are 6.9 x 1(T 6 (7.7 x 10- 5 ), 2.5 x 1(T 5 (2.6 x 10" 4 ), and 6.0 x 1(T 6 (3.2 x KT 4 ), and that 
of D° — > / + /~7 (I = e, fx) and D° — > z/z/7 are 6.3 x 10 -11 and 2.7 x 1(T 16 , respectively. 



1 Introduction 



In the standard model, flavor changing neutral current (FCNC) decays of D mesons 
are predicted to be very rare. Such processes are forbidden at the tree level but can occur 
at higher loop level. On the other hand, these rare decays make a good chance to probe 
for new physics M. Recently, the new experimental limits with the range 10~ 4 to 10~ 5 
for FCNC decays of the charm mesons were obtained at Fermilab E791, E771 and CLEO 

ni 1,1. 

As we known, the two-body leptonic decays of D mesons of D — > ijj (l it j = e, //, v) 
are helicity suppressed. Explicitly, in the standard model, the decay branching ratios of 
D ( S ) e+z/ e, r+u T andL>° -> l+l~ are l.OxlO -8 (9.5xl(T 8 ), 4.7xl(T 4 (9.5xl(T 3 ), 

1.2 x 10~ 3 (4.2 x 10 -2 ) and O(10" 19 ), respectively, while that of D° —> vv is zero. It is 
clear that the rates for the modes of D — > 1^ (Zj = e, /i, v) are too small to be measured 
due to the helicity suppression as well as Glashow-Iliopoulos-Maiani (GIM) mechanism 
effect 0. 

To avoid the helicity suppression, we will consider three body decays such as D + — > 
/ + z/;7 and D° — > Z + Z~7 with a photon radiated from charged particles. The ampli- 
tudes of the decays can be divided into the "internal-bremsstrahlung" (IB) parts and 
the "structure-dependent" (SD) parts. Since the photon is emission from the external 
charged leptons, IB parts are still helicity suppressed for the light charged lepton modes, 
while in SD one of the photon is emitted from intermediate states that they are free of the 
helicity suppression. Therefore, the decay rates of D + — ► / + z/;7 and D° — > / + /~7 (I = e, fi 
) might have an enhancement with respect to the purely leptonic modes of D + — > l + ui 
and D° — > , respectively. 

It has been pointed out that the heavy flavor decays can be analyzed with the help 
of heavy quark effective theory. For charmed decays, this is questioned since the charm 
quark mass is not much larger than the QCD scale. Therefore, we will use the light 
front quark model || [7|, |8], |9], 10, |ll[ to evaluate the hadronic matrix elements in D —>■ 7 



transitions. These decay modes have been studied in various models [12, 13, 14]. It is 
known that as the recoil momentum increases, we have to start considering relativistic 
effects seriously. The light front quark model |7| is the widely accepted relativistic quark 
model in which a consistent and relativistic treatment of quark spins and the center- 
of-mass motion can be carried out. In this work, we calculate form factors in D — > 7 
directly at time-like momentum transfers by using the relativistic light-front hadronic wave 
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function. The parameters in the wave function are determined from Refs. J7|, Within 
the framework of light-front quark model for decay processes, one may calculate hadronic 
matrix elements in the frame where the transfer momentum is purely longitudinal, i.e, 
p± = and p 2 = p + p~ , which covers the entire range of the momentum transfers. We will 
give their the momentum dependence in whole kinematics region of < p 2 < p^ax- 

The paper is organized as follows. In Sec. 2, we study the form factors of D — > 
7 transitions within the light front framework. We present the decay amplitudes and 
calculate the decay branching ratios of Dt\ — > / + z^7, D° — > vv^i and D° — > in 
Sec. 3, 4 and 5, respectively. In Sec. 6, the conclusion is given. 

2 Form factors of D —> 7 transitions 



To evaluate the D u — > 7 transition , one needs to consider the following currents: 

J M = U7"(l-75)c 

= uia^p u (l + l5 )c, (1) 

where p is the momentum transform, D u = uc and u stands for the light quark such as 
a d-quark or u-quark or s-quark so that = D + , D u = D° and D s = Df, respectively. 
Our main task is to calculate the hadronic matrix elements for D — > 7 transitions, which 
can be parametrized as: 

(P • <? 7 )< - ( e * ■ P)Qtt] > 
(l(q 7 )\u-f^c\D(p + qi )) = iej^-e^e^pPq!^, (2) 

('j(q 7 )\uia flu p u 'j 5 c\D(p + g 7 )) = ieF TA e u (p ■ q^g^ - p v q~ m ), 
('~l(q~ i )\uio lJLP p v c\D(p + g 7 )) = eF TV e^i3pxe p qPp x , (3) 

where the e M is the photon polarization vector, g 7 and p + q 1 are photon and D meson 
four momenta and Fa, Fy, Fta and Ftv are the form factors of axial-vector, vector, 
axial-tensor and tensor, respectively. 

In the light front approach, the physically accessible kinematics region is < p 2 < 
Pmax = due to the time-like momentum transfers. The D u meson bound state in light- 
front quark model consists of an anti-quark u and a quark c with the total momentum 



(l(q~f)\uivl5c\D(p + q^)) = -e— 
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(p + Q-r)- ^ can be expressed by: 

\D(p + q 7 )> = £ J[dk l ]{dk 2 }2{2'Kf5\p + q ry -k 1 -k2) 



A1A2 ' 



x $£> lA2 (x, ki_)b+{h, Xi)dt(h, A 2 )|0 > , (4) 



where ku2) is the on- mass shell light front momentum of the internal quark u(c), the light 
front relative momentum variables (x, k±) are defined by 

k+ = x(p + g 7 ) + , k 1± = x(p + g 7 ) ± + k± , (5) 

and 



2k + k + \ 3 

K lX2 (x,k ± )= ( 12 -3] 11(^1)7^,^)^*1), (6) 

.-Mq — \ m u — m b) 



with k±) being the momentum distribution amplitude. The amplitude of <f)(x, k±) can 
be solved in principles by the light-front QCD bound state equation 0, §, ||. However, 
we use the Gaussian type wave function: 



(l , W = i NM_^exp ], (7 ) 



where 



r ,, n dk + dk± „ T i / 7r " 

■*J = W iV = 4 U, 



For the gauged photon state, one has: 



| 7 (<? 7 ) > = iV'ja+(g 7 , \) + J2f [dk l }[dk 2 }2{2iif5\q 1 -k x - k 2 ) 
I AiA 2 J 



where 



x $^ A2A (g 7 , fci, fe 2 )6+(A;i, A x )4(A:i, A 2 ) |0 > , (9) 



7_L • fci - mi 



7-1 • e± ? XAi , (10) 
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with 

ED = ^_kl+ml_kl + ml_ 
q+ kj kj 

To calculate the matrix elements of the Eqs. @ and (f|), we choose a frame with the 
transverse momentum p± = 0. The p 2 = p + p~ > covers the entire range of the mo- 
mentum transfers. By considering the "good" component fi = +, the hadronic matrix 
elements of Eqs. @ and @ can be rewritten as: 

p 

<l{q 1 )\uXlsC+\D(p + q 1 )> = ~ ie 7^( e l -%±)P + i 

< 1 (q 1 )\uXc + \D(p + q 1 )> = e ^-e ij e*q ljP + , (12) 

and 

< 7(9 7 )I( M +7°75C- - utl°l 5 c + )\D(p + q 7 ) > = eF TA (e* • q 1± ) , 

< 1 {q 1 )\{u% 1 Q c^-utl Q c + )\D{p + q 1 )> = - ie F TV ^e\q^ , (13) 

where w + (c + ) and w_(c_) are the light-front up and down components of the quark fields, 
and they have the two-component forms: 

d+ = ( o 1 > ( 14 ) 



and 



(iaj_ ■ d± + f3m q )q + 



id+ 
( ° 



(15) 



respectively. 

In Eq. (15), Xq is a two-component spinor field and a is the Pauli matrix. The form 
factors Fy,A and Ftv,ta m Eqs. (|l^) and flHf) are then found to be: 

f2 m b -Bkl<d 2 m u -Ak 2 ± Q \ 
X \3 m 2 c + kl 3 ml + ki J ' 1 j 



2 m c + (1 - x')(m c - m u )k]_Q 2 m u - x' (m c - m u ) k 2 ± Q \ 

3 m 2 + ^ + 3 m 2 + fc 2 | ' 1 j 
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2(2tt) j 



[2A 1 + A 2 A;ie 2B l + B 2 k\Q 
x <; — h 



3 m 2 + k 2 ± 3 ml + k 2 ± 



F TV (P 2 ) = ~ J 



dx d 2 k 



2(2 



7T 



w $(x',k 2 ± 



( 2 d + C 2 k 2 & 2D 1 + D 2 k 2 & 

x < - — - — — h -- 



3 m 2 + /c 2 



3 + A;^ 



(1 — 2x)x'(m c — m u ) — 2xm u , 

x' + x - 2x'x) [x'(x - 1) - x(2x - 1)] k]_ 



2(1 - x)x'm c + (1 - 2x)(l - x')m u 
2 



xx l2 {l -x')(l-x) 
+x [(x — x') + 2x'x(l — x)] m 2 + 2x 2 (l — x') 2 m c m u j , 

+2x(l — x')m u m c — x{l — 2x){l — x') 2 m 2 + x'(l + x'x — 2xx' 2 )m; 

— -J (x' + x- 2x'x)(l -2x + 2x 2 - x'x) A; 2 

x'x(l -x)(l -x') 2 \ V A ; x 

+2xx'{l - x) (1 - x)m u m c + (1 - a/) (x' + x - 4x'x + 2x'x 2 )m 2 j 

2(x - x') 



— <^ - (x' + x- 2x'x)(l - 2x)/c 2 
x'x(l-x)(l-x') 2 \ V 

— "(1 - 2x)x' 2 (l - x)m 2 + (1 - x')(x' + x - 3x'x - 2x 2 (l - x'))m\ 
2 

xx' 2 (l-x')(l-x) 



(x — x' + xx')(x' + x — 2x'x)k 2 



— 2x 2 (l — x') 2 m c m u + x' (x — x') — 2(1 — x')x 2 

2(X — X ) J . n /„\,2 i / 



xx' 2 (l -x')(l -x) 



(x' + x — 2x'x)k ± + x'(l — 2x + xx')m ( 



—m q x(l — x') — 2x(l — x')m u m c > , 



(1 - x)(l - 2x + x')(x' + x - 2x'x)A; 2 



x(l — x)x'(l — x') 2 k 
- (1 - x')(x' + x - 2x 2 - 2x'x + 2x 2 x')m 2 + 2x' 2 (l - x) 2 m u m c 



D 2 = 

x 



'(l — x)x(i — x'y y 



+ (1 — x')(x' — x — x'x)m^j , 

1/2 



$(x,fci) = iV ( t ^ \ -r^exp 



Mq — (m u — m c ) 2 ) \ dx \ 2uj 2 d 



_ 1 d®(x, k\ 



$(x, kj_) dk\ 



x = x(l-j^J, k=(k ± ,k z )- (20) 

To illustrate numerical results, we input m u = 0.3, m c = 1.4, and to = 0.5 in GeV. 
The form factors of Fy,A and Ftv,ta are shown in Figures 1 and 2, respectively. 



3 Ifc-Z + H7 



3.1 Matrix elements and kinematics 



In this work, we investigate the decays of — > (q = d,s). The effective 

Hamiltonian for c — > g at the quark level in the standard model has been obtained by 
using similar results for b — > g Zz// 

#e//(c -> <?^) = -^Kg<hv(l - 75)c^/7/,(l ~ 7s)i • (21) 

From Eq. due to the emissions of real photons that the decay amplitude of D+ — > 

/ + z//7 can be written as 

M = M IB + M SD , (22) 



M /B = l ^V cq f D m l elD^ (23) 

M 5D = -i^LVcqeaLpH * , (24) 

where Mjb and M$d represent the amplitudes from the IB and SD contributions, respec- 
tively, and 

D* = u(p u I+75 ^^)v(p llSl ), (25) 

p • q 2pi ■ q 

L/3 = u{p v )^fp{l-^)v{p h si), (26) 

RaP = ^(-9 a(3 p-q + P a q l3 )+ie-^e al3 ^q,p u , (27) 



with /d, being the D( s ) decay constants, e M and si the four vectors of the photon po- 
larization and charged lepton, and p, pi, p v , and q the four momenta of Df s \, I, v and 7, 
respectively. The physics allowed region of Df-, — > / + z^7 is 



mf <p 2 < M 2 D 



(2f 



To describe the kinematics of Dfs — > / + z^7, one needs two variables, for which we define 
x-y = 2E 1 /M D(g) and y = 2Ei/M D{s) . We write the physics region for x 7 and y as: 



where 



1 - x~ + 

1 — x^ 



< x 7 < 1 — r; 
< y < 1 + n , 



mf 
W D 



7.1 x 10~ 8 

3.2 x 10- 3 
0.9 



for I — e 
for I = p. 
for I = t 



(29) 
(30) 



(31) 



for D + — > / + ^7 and 



?7 



6.4 x 10~ 8 for I = e 
2.8 x 10~ 3 for l = fi 



0.8 



for I = t 



(32) 



for — > l + vi r ). 



3.2 Decay rates 



In the (q = d or s) rest frame, the partial decay rate is found by |TJ| 



1 1 



(2vr) 3 8M 



M I 2 d£Ld£, 



Using x-y and y, from Eq. ([J3|) we obtain the differential decay rate as 



and 



Mo ...„ M D ^ (1 _ A)A 



dx^dy 



256tt 3 



IMI 



256tt 3 



(33) 



(34) 



A = A /B (x 7 , y) + A SD (x 7 , y) + A /iV (x 7 , y) , 



4m 2 1 /. 



D a 



x 2 + 2(1 - n) ( 1 - x 7 - ^ 



(35) 
(36) 



7 



\F V + F A \ 2 



I- XV 7 A 



+ \F v -F A \\y-X) 



(37) 



A IN {x 1 ,y) 



-4M Dq mi 
-Re[f D (F v -F A 



Re[f Dq (F v + F A )*\[\ 



_ n 



A 



where A = (x 7 + y — 1 — r)/x 1 . In Figures 3 and 4, we show the differential branching 
ratios of -D + — > /i + z^t7 and £)+ — > p + i/^7 as functions of the photon energies, where we 



0.22, \V a 



0.974, 



have used rrid = 0.3 GeV, m s = 0.4 GeV, m c = 1.5 GeV, \V c d 
uj = 0.5 GeV, r D + ~ 0.467 ps, r D+ ~ 1.05 ps and / D = = 230 MeV p|, 0. We g 
the integrated branching ratios of D + and — > l + vi r y (I = e, p, r) in Tables 1 and 2. 

Table 1: Integrated branching ratios for the radiative leptonic D + decays 



Integrated Decay 










Ref. 


Branching Ratios 


IB 


SD 


IN 


Sum 





10 6 B(D + -»■ e + z/ e7 ) 


1.0 x 10~ 3 


6.9 


6.2 x 10~ 5 


6.9 


82 


10«B(D+ -> p+^7) 


17.2 


6.8 


6.2 x 10- 1 


24.6 




10 6 £(£>+ -> t + z/ t7 ) 












(Cut 5 = 0.01) 


5.6 


1.3 x 10~ 6 


3.7 x 10- 1 


6.0 





Table 2: Integrated branching ratios for the radiative leptonic Df decays 



Integrated Decay 
Branching Ratios 


IB 


SD 


IN 


Sum 


Ref. 

14| 


Ref. 

[13| 


10 5 5(£>+ -> e + u e j) 


1.1 x 10~ 3 


7.7 


6.8 x 10~ 5 


7.7 


90 


17 


10 5 B(Df - n+vft) 


17.8 


7.5 


7.1 x lO^ 1 


26.0 




17 


10 b B{Df -> r+z/ r7 ) 
(Cut 5 = 0.01) 


30.7 


1.8 x 10~ 3 


1.7 


32.4 







We now compare our results with Refs. [13], 0|. As shown in Tables 1 and 2, our 



result are much smaller than that in Ref. [14| but close to the values in Ref. [ 13| . We 
note that for the p and r channels the IB contributions are dominant since D-mesons are 
not very heavy. 
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D° z/z/7 



4.1 Matrix elements and kinematics 

For the processes of D° — > v\Vi^ (I = e,fi,r), at the quark level, they arise from the box 
and Z-penguin diagrams that contribute to c -> uv\V\ with the photon emitting from 
the charged particles in the diagrams. However, when the photon line is attached to the 
internal charge lines such as the W boson and quark lines, there is a suppression factor 
of m^/M^r in the Wilson coefficient in comparing with those in c — > uv\Vx. Thus, we need 
only to consider the diagrams with the photon from the external quarks. The effective 
Hamiltonian for c — > uv\V\ is given by 

G ex 

H eff(D° -> Viva) = — g- — r^— Yl V ui V*D(xi)u-f^l - 75) a>a v {l - 75)14 , (39) 

V2 27rsm V w i=d b 



where 



D{xi) = \ 



xf-Xj + Z 3(^-2) 

Xi—1 (1 - XiY 



(40) 



and Xi = m^/M^r with rrii (i = d, s, b) being the current quark masses. We note that in 
Eqs. (^|) and (fiO"D , only the leading contributions have been included and the additional 



a, correction to the result, which are small [18 



4.2 Decay amplitude 

From the effective Hamiltonian for c — > uv\V\ of Eq. (|39"D and the form factors defined in 
Eq. fljjD , we can write the amplitude of D° — > vivf) as 

M = -*H L a a 2« E V m V*D(^KH^u(p p ) ltl (l- l5 )v(Pv), (41) 

with 

F F 
Hf,u = jt-(-P- Q~< fa + PvQjv) + itnvaf3^q"p P • (42) 

1V1 D Md 

where Md is the mass of D° and the form factors are given by Eqs. flTE| ) and (|T7|). 

Similar to the decays discussed in the previous subsection, we also define x 1 = 2E 1 /Md 
and y = 2E d /Md in the D-meson rest frame in order to re-scale the energies of the 
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photon and anti-neutrino. By integrating the variable y in the phase space, we obtain the 
differential decay rate of D° — > z/z/7 as 

where we have included the three generations of neutrinos. 

Using m c = 1.4 GeT^, m u = 300 MeV, and uo = 0.5, the differential decay branching 
ratio dB(D° — > vurf/dx^ as a function of x 7 = 2E 1 /Md is shown in Figure 5 and the 
integrated decay branching ratio is 

B(D°^uuj) = 2.7 xl0~ 16 . (44) 

It is easy to see that in the standard model the branching ratio of D° — * z/z/7 is too 
small to be measured experimentally. However, a large branching ratio may be induced 
by new physics such as the leptoquark model |Lj| . 

5 D° /+/- 7 



5.1 Matrix elements and kinematics 



The contribution to the process of D° — > l + l~ r y (I = e or //) arises from the effective 
Hamiltonian that induces the pure leptonic mode of -D — > Z + Z~. The short distance con- 
tribution for c —>■ u comes from the W-box, Z-boson and photon penguin diagrams, 
and the effective Hamiltonian is given by 



2V2irsm z e w i=d b { 



Ff sin 2 6 



w 



p- 



m c uia flI/ p"P R cl>y fl l 



(45) 



where Pl,r =(1 T7s)/2, and p is the momentum transfer which is equal to the momentum 
of the lepton pair. The Wilson coefficients F/, F 2 , F 3 can be found in Refs. 0, ^] and 
one has 



F/ = + C? - sin 2 MQ 5 + 2C? 



F 



-sin 2 %(C? + 2C 2 



+ 



-Q 
1 1 



1 1 



+ 



Axi-l 4(^-1) 
9 1 3 



2 (x* - l) 3 



3 Xa In 



2 (x< - l) 4 



2x»-l 4(x i -l) 2 2(^-1) 



3 x 3 lnx 7 - 



^2 r\ 



2 



(46) 
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where Xj = rn^/M^ with m, (i = d, s, b) being the current quark masses. The coefficients 
of Cj, Cf and Cf are given by 



C} 



c 



1 



Xj ill 



— 1 (Xj - l) 2 
Xj 3 1 3 2;XL- Xt 



(Xi - I) 



11 13 1 

+ 



4 8 Xi - 1 





7(6, Xi) , 

;\nxi + 7(f,Xi) 
1 1 



,1 ; 



+ 



YlXi-\ 12(xi-l) 2 2(x;-l) 3 . 

1 1 



2 


1 




+ 


(1 




1 




5 




1 




3 Xi 




1 








I) 2 


6 


(Xi 




1)3 


7 


1 




+ 


13 




1 




1 




1 




3 Xj 




1 


12 


yX% 




I) 2 


2 


yXi 




1)3 


1 


1 






35 




1 




5 




1 




6xi 




1 




12 


[Xi 




I) 2 


6 


[Xi 




1)3 



In Xj 



+ 



1 1 



Xilnxi - 2y(£,Xi) , (47) 



and 



7(C,Xi 



1 /3 1 



1 1 



£xj — 1 \4xi — 1 8£xj — 1 

1 V+i 



I X,; In X,; 



1 1 1 



• £ - 1 



(48) 



where £ is a gauge dependent parameter. The decay amplitude for D° — > l + l~y can be 
written as 



M(D° -> /+r 7 ) = A4 /B + M 5 d 



(49) 



where .M/s and -Ms_d represent the IB and S.D contributions, respectively. For the IB 
part, the amplitude is clearly proportional to the lepton mass mi and it can be written 
as: 

GpQ. 

2^2tt sin 2 8 W A , 



M 



IB 



-ie- 



E v ul v;j Dmi (F^ - fI 



D 



-Hi 



,(50) 



2pi-g 7 2p 2 -q~ / ' 

where Pd, Pi, P2 and g 7 , are the momenta of D, l + , l~ and photon respectively. The 
amplitude for the SD part, where a photon is from the initial quark line, is given by 

Gpoi 

■ a 

i=d,s,b 



M 



SD 



, T . ,„ E VmV* \Ffa(q y )\u^P L c\D(p + qJ)l'f(l-'y 5 )l 



+ F?(y(q,)\u lfl P L c\D(p + g 7 ))/ 7 M (l + 75)/ 

- Z ^'f ^ {l{q,)\uia»uP u P R c\D{p + q^l 



(51) 



which shows that to find the decay rate, one has to evaluate the hadronic matrix elements 
in Eq. (KB. 
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Using Eqs. (2) and (3), we rewrite Eqs. Q5TD as: 

Gpa 



M 



SD 



2 v / 2vrsin 2 ^ 



W i=d,s,b 



+ I 



e!(p-g 7 )-(e*-p)g 7/1 Bl^l + Dh^l 



where the factors of A-D are 



A = 


(Fl + F 2 ) 
M D 


-F VA (p 2 ) - 


2F^F TA (p 2 


£ = 


(Fl + FT, 
M D 


-F VV (P 2 ) - 


2F?^FMP 2 


c = 


(F 2 - Fl) 
M D 


-F VA (P 2 ), 




D = 


(fi - fi: 


-Fvv(p 2 ), 





respectively. 

5.2 Decay rates 

The partial decay width for D° — > l + l~ r y in the D° rest frame is found to be: 
1 dq dpi dp 2 



dT 



.\M\ 2 (2ny5(P- Pl -p 2 -q,)- 



2M D r ^ a,/ (2vr) 3 2£ 7 (2vr) 3 2£i (27r) 3 2£ 2 ' 

where the square of the matrix element is given by 

\M\ 2 = \M IB \ 2 + \M SD \ 2 + 2Re(M IB M* SD ) . 

To describe the kinematics of D° — > / + /~7, we define two variables of x 7 = 2P B 
and y = 2Pp ■ pi/M^. We write the transfer momentum p 2 in term of x 7 as: 

p 2 = M 2 D {l-x 1 ). 

The double differential decay rate is then found to be: 

d 2 T l M, 



dXjdX 



256tt 3 



L> \M\ 2 = Cp(x„X) } 



where 



C = a E \^^LfG 2 F Ml 



i=d,s,t 



and 



p(x 7 ,A) = P/b(x 7 , A) + Psd(x^, A) + p/ A r(x 7 , A), 
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with 



PSD 



Mix* 



A(l-A) 



2r,(^ 



A 2 (1-A) S 



.1/ 



(|A| 2 + (1 - x 7 + 2n) - 2(1 - x 7 )(A - A ) 



+(|Cf + |D| 2 ) (1 - x 7 - 2n) - 2(1 - z 7 )(A - A 2 ) 



+2Re(B*C + A*D)(l - a; 7 )(2A - 1) j, 
P/iV = /D(^ 2 -^)ri/?e(A)— ^— + i?e(D)-^ (1 " 2A) 



(60) 



A(l-A) v ' A(l-A) J ' 

Here A = (x 7 + y — l)/x 7 and r; = m 2 /M% and the physical regions for x 7 and A are given 
by: 



o < Xr. < l - An 



1 - -.1 

2 2 V 



4n 



1 — a;- 



< \<- + -Jl- 
-2 2V I 



4n 



(61) 



In Figure 6 we present the differential decay rate of .D — > /i + /i~7 as functions of 
x 7 . From the figures we see that the contributions from soft photons, i.e the IB ones, 
corresponding to small x 7 region, are infrared divergence. To obtain the decay width of 
D° — > / + / _ 7, a cut on the photon energy is needed. The integrated branching ratios are 
summarized in Table 3. Here, we have used the cut value of 5 = 0.01 and m c = 1.5GeV, 
m s = 0.4 GeV and = 0.3 GeV in the calculations of the form factors. We have also 



used \V cd V* d \ = \V CS V* S \ = 0.2145, \V ch V* b \ = 1.6 x HT 4 , t(D°) = 0.415 x 10" 12 s @ to 
evaluate the numerical results. 



Table 3: Integrated branching ratios for the radiative leptonic D° decays 



Integrated Decay 
Branching Ratios 


IB 


SD 


IN 


Sum 


10 U B(D° -> e+e- 7 ) 


1.3 x 10~ 9 


6.3 


4.6 x 10" 4 


6.3 


10 U B{D° -> n+yri) 


1.9 x 10~ 5 


6.3 


2.2 x 10~ 3 


6.3 



We see that the decay branching ratios of the radiative leptonic decays of D° — > / + /~7 
are about eight orders of magnitude larger than that of the purely leptonic modes of 
D° — > for I = e and //, respectively. 
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6 Conclusions 



We have calculated the form factors of the D — > 7 within the light-front model and 
used these form factors to evaluate the branching ratios of the leptonic .D-meson radiative 
decays. We have found that, in the standard model, the decay branching ratios of 
e+z/ e 7, Djj -> // + iyy and -> r + z/ r7 are6.9xl(r 6 (7.7xl0~ 5 ), 2.5 xKT 5 (2.6 xlO" 4 ), 
and6.0xl0" 6 (3.2 xlO -4 ), and that of £>° -> /+/"7 (/ = e,//) and D° -> ^7 are6.3xl0- n 
and 2.7 x 10~ 16 , respectively. Comparing with the purely leptonic decays of D° — > , 
we have found that the rates of the corresponding radiative modes are several orders 
of magnitude larger. We conclude that some of the radiative leptonic decays could be 
measured in the BTeV/CO experiment at Fermilab [[Hp. During Tevatron Run II, which 
could reconstruct about 10 9 charm mesons decays and increase the upper limit of about 
three orders. 
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Figure Captions 



Figure 1: The values of the form factors F v (solid curve) and F A (dashed curve) as 
functions of the momentum transfer p 2 for — > 7. 

Figure 2: The values of the form factors F T v (solid curve) and F T a (dashed curve) as 
functions of the momentum transfer p 2 for Df — > 7. 

Figure 3: The differential decay branching ratio dB(D + — > n + v^) j 'dx 7 as a function 
of x y = 2E 1 /M D . 

Figure 4: The differential decay branching ratio dB(D+ — > ii + u^) / dx 1 as a function 
of x 7 = 2E 1 /M Ds . 

Figure 5: The differential decay branching ratio dB(D° — * vv^jdx^ as a function of 
x 7 = 2E 1 /M D . 

Figure 6: The differential decay branching ratio dB(D° — > fi + fi~^y)/dx y as a function 
of x 7 = 2E 1 /M D . 
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Figure 1: The values of the form factors Fa (solid curve) and Fy (dashed curve) as 
functions of the momentum transfer p 2 for — > 7. 
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p 2 (GeV 2 ) 

Figure 2: The values of the form factors Fta (solid curve) and Ftv (dashed curve) as 
functions of the momentum transfer p 2 for — > 7. 
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